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Abstract 



Using a recently developed perturbation theory for uasinormal modes 
(QNM's), we evaluate the shifts in the real and imaginary parts of the QNM 
frequencies due to a quasi-static perturbation of the black hole spacetime. 
We show the perturbed QNM spectrum of a black hole can have interesting 
features using a simple model based on the scalar wave equation. 
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I. INTRODUCTION 



The observational consequences of black holes interacting with their astrophysical envi- 
ronments have been a subject of much interest for the last 30 years. Within the next few 
years, it is expected that the new generation of gravitational wave observatories (LIGO, 
VIRGO) will be able to detect gravitational waves emitted by black holes excited by 
matter, or even other black holes, falling into them. It has long been known that the 
gravitational waves emitted in such a process will carry a signature associated with the 
well-defined quasinormal mode (QNM) frequencies of the black hole 0, and will, among 
other things, provide confirmation of the existence of black holes. Numerical simulations ]3[] 
suggest that in some cases the QNM ringing may even dominate the signal. 

QNMs of black holes have been extensively studied with the black hole perturbation 
theory ||. If a black hole settles down in an otherwise empty and asymptotically flat 
spacetime at the end-point of dynamical evolution, it will be a Kerr black hole (Schwarzschild 
hole in the case of zero angular momentum) |J. Weak (linearized) gravitational waves 
propagating on the Kerr or Schwarzschild background can be described by the Klein-Gordon 
equation ||: 

—kg 

d?-d? + V(x)]$(x,t) = 0, (1.1) 

where x is a radial coordinate, $ is the radial part of a combination of the linearized changes 
in the metric functions representing the gravitational wave, and the outgoing wave boundary 
condition is appropriate for waves escaping to infinity. The potential V(x) describes the 
scattering of the gravitational waves by the background geometry. For example, in the case 
of a Schwarzschild hole of mass M, V is the Regge- Wheeler potential MM ■ — schw 



V(x)=(l-^) 



/(Z + l) 2 ,2M 



;i-2) 



for each angular momentum sector I, where x = r + 2Mln (r/2M — 1), s is the spin of the 
field (s = 2 for gravitational waves), and r is the circumferential radius. 

A single-frequency solution [$ oc exp(— iut)] with the outgoing wave boundary condition 
is a QNM, with Im uj < 0. The QNM spectra of Kerr and Schwarzschild black holes have 
been extensively studied theoretically Q and numerically and the positions of the 

QNM frequencies in the cu-plane are known in detail ]8|I|. 

For each given / in (|1.2j ), the QNMs extend downwards in a string in the tu-plane, with 
Re uj nearly constant and Im uj nearly uniformly spaced || (Fig. [I|). The known pattern 
of frequencies provides a template against which one can try to determine the nature of 
the source. For an isolated black hole, the no-hair theorem || implies that the spectrum 
is described by only two parameters, the mass M and the angular momentum J of the 
hole. However, the black holes that are likely to be observed will not be isolated, but will 
likely be situated at the centers of galaxies, or will be surrounded by massive accretion 
disks. Therefore the observed spectra should not be matched against those of a pure Kerr or 
Schwarzschild black hole, but to one perturbed by its surrounding — a dirty black hole. We 
should immediately caution the reader that while gravitational waves from black holes are 
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expected to be detected within the next few years, a determination of the QNM spectrum 
with the frequencies of a few modes included might not be an easy task. Indeed to what 
extent the gravitational radiation from realistic black hole events would be dominated the 
QNM spectrum is still a matter of much controversy. However, in as much as the goal of 
the gravitational wave observatories is to obtain astrophysical information of our universe 
(the "O" in "LIGO"), there is no doubt that we will eventually have to face this problem of 
the QNM spectra of dirty black holes. 

We note that two kinds of perturbations are involved here. In the standard black 
hole perturbation theory ||, ( |1.1| ) is obtained by linearizing the metric about the Kerr 
or Schwarzschild background, and the time-independent eigenvalue problem (with the out- 
going wave boundary condition) determines the QNM spectra of isolated holes. The second 
type of perturbations are the perturbations that change the background on which the wave 
propagates, e.g., by the presence of an accretion disk near the black hole. The perturbation 
of the background can often be regarded as quasi-static, and hence separable from that of 
the gravitational wave perturbation by the time scales involved (in a suitable gauge choice). 
In this paper we focus on time-independent perturbation of the background described by 
( |1.1| ) with a potential V(x) = Vo(x) + /iVi(x), \fi\ <C 1 (a model problem is given in Sect. [TV| 
below). We are led to study the QNM frequencies of the following eigenvalue problem in 
powers of \i: 

-kgl 

- (j>"{x) + [V (x) + fiV^x)] <f> = u 2 <t>. (1.3) 

(|TT3|) is appropriate for considering the perturbed Klein-Gordon wave equation describing 
the propagation of scalar waves in a gravitationally perturbed (dirty) black hole spacetime 
(as will be shown below). In this paper we show how the disturbed QNM spectrum can 
be determined for such a system. This represents the first step towards determining the 
disturbed gravitational wave QNM spectrum of a dirty hole, the real case of physical interest. 
For a realistic black hole perturbed by an external matter fluid source, the gravitational wave 
QNM spectrum also involves the fluid modes. Just as in the case of the perturbation of fluid 
star, we would expect two types of perturbations, with one strongly involving the fluid 
motion (e.g., the polar f,g,p modes), while the other only weakly involving the fluid (e.g., 
the axial modes). We expect the scalar perturbation studied in this paper to be more easily 
generalizable to the latter kind of modes. A more complicated set of equations would have 
to be used to describe the former type of perturbation involving the "fluid" modes in the 
shell of matter outside the hole. While the calculation in this paper may not be easily 
generalizable to these more complicated situations involving gravitational waves coupled to 
matter, we note that a full perturbation treatment of the gravitational wave case will not 
be possible without a thorough understanding of the behavior of the Klein-Gordon wave 
equation with a perturbed potential, namely the system studied in this paper. 

While the perturbed Klein-Gordon wave equation ( |1.3| ) is superficially similar to stan- 
dard textbook problems, e.g., the usual Rayleigh-Schrodinger perturbation theory (RSPT), 
we note that the perturbation problem encountered here is fundamentally different: the 
outgoing wave condition renders the system physically nonconservative (energy escapes to 
infinity) and the associated operator [— d 2 /dx 2 + V(x)} non-hermitian; hermiticity underpins 
the usual RSPT. 
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The difficulty coming from the non-hermiticity can be seen in several guises if one tries 
naively to transcribe the usual formulas. Unlike the hermitian case, now the unperturbed 



eigenstates do not in general form a complete set for expansions [rO , IT| , at least not in 
the case of black holes. The usual RSPT formula in terms of a sum over intermediate 
eigenstates is therefore inapplicable. Even the first-order shift, which does not involve a 
sum over intermediate states, cannot be given by the usual formula {^>q\^iV\\4>o) / (0o|0o), in 
obvious notation — the usual inner product leads to (0o|</>o) = /^oo^^o^o = 00 since a 
QNM wavefunction extends over all space (and indeed grows exponentially at infinity). 

So far, the perturbation of black hole QNMs has attracted little attention, partly be- 
cause a perturbative formalism for the QNMs of an open system, as opposed to the normal 
modes (NMs) of a conservative system, has not hitherto been available. In this paper we 
develop such a formalism, which then opens the way to extracting information about the 
astrophysical environment of the black holes from the observed signal, beyond the mass and 
the angular momentum of the hole. 



This paper is a follow-up of fi~2f , which outlined some of the results derived in this 



paper. In Sect. [IT] of this paper, we develop a formulation for the perturbation of QNM 
systems. As a first step in this direction, we limit ourselves to the scalar wave case, in which 
the evolution is described by a single Klein-Gordon equation with a perturbed potential. 
The shifts in both the real and imaginary parts of the QNM frequencies u are obtained in 
quadratures in terms of fiVi(x), in principle to arbitrary order in [i. Given the precision of 
the observational data that is possible in the near future (indeed at this point it is not clear 
how many QNM's one can extracted from the waveforms of black hole events, given the 
S/N of even the advanced phase LIGO), the emphasis is on the first-order shift. The shifts 
when expressed in terms of a generalized inner product take a form similar to that in RSPT. 
The perturbative results for a Schwarzschild black hole are derived in Sect. [TTT[ We show 
that a function H(x) can be defined which depends only on the original unperturbed system. 
We investigated and presented in detail the properties of this function H for the black hole 
case. This function controls the phase and magnitude of the first-order shift of the spectrum 
for any given perturbation (not just for the model problem in this paper), hence providing 
insight to the properties of the black hole spectra in general. Sect. [TV] illustrates these 
results with a model problem where the perturbation is due to a spherical shell of matter 
located at a fixed radius, and we study scalar wave propagating in this background. The 
frequency shifts are obtained using the perturbation formula, and compared to numerical 
results. We show even in this simple case the spectra for shells located at different radii 
contain very interesting features. These features can be understood by the perturbation 
formula, demonstrating the power of the perturbation formula in providing understanding 
of the perturbed spectra. 



II. FORMULATION 



Waves defined by differential equations such as ( |1.3|) or the optical analogs described 
by the wave equation are open systems if the outgoing wave boundary condition is 
imposed; as such they are physically nonconservative and mathematically non-hermitian. 
The usual tools of mathematical physics based on the hermiticity of the defining operator 
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would not, in general, be expected to apply. In attempting to develop generalizations of 
the familiar formalisms for conservative systems, it is important to recognize that these 
open systems fall into two broad classes. In the first, the potential V(x) has discontinuities 
at x = 01,02, and vanishes at infinity suitably fast ("no tail"); when these conditions are 
satisfied, the QNMs form a complete set on the interval [01,02], and the usual formalism 
can be carried over with minimal changes |lO| Jl^| . This is in itself somewhat surprising, in 
that the dynamics are controlled entirely by the resonances without the need to add any 
"background" . However, when the discontinuity or "no-tail" conditions are not satisfied - 
and this is the situation of interest in this paper — the QNMs are in general not complete, 
and the dynamics is not completely controlled by the QNMs. The QNMs are manifested 
in the complex frequency plane as poles of the Green's function; now, there is in addition 
a cut along the negative Im uj axis fI5|| . The cut arises from scattering by the asymptotic 
part of the potential (in the sense that if V(x) is truncated at any finite distance, the cut 
disappears). This is of course the case for the radial problem of a black hole, where V(x) has 
no discontinuity and goes asymptotically as a centrifugal barrier plus ~ logx/x 3 ; this tail 
generates the cut in the complex frequency plane, extending all the way to Im oj — ► 0, and 
consequently leads to the long time power- law behavior ~ t~ a of the dynamics [|15[|16| . The 
point to be stressed here is that since the QNMs are not complete, a perturbation formalism 
based on a sum over intermediate states will not be appropriate. 

In this paper we need to rely on a different approach, which does not require a complete 
set of QNMs, by generalizing the logarithmic perturbation theory (LPT) |T7|-f2~4|j , which 
focuses on the logarithmic derivative f(x) = <p' (x) / 4>(x) . Three features are relevant in 
the present context. First, the numerical determination of QNMs is known to be difficult, 
essentially because the exponentially growing solutions are numerically sensitive; this makes 
it important (indeed more important than would be the case for NMs) to develop analytic 
techniques, as well as semi-analytic techniques such as perturbation theory. Second, LPT 
does not require a complete set of unperturbed states for expansion, and is therefore well 
suited to situations where QNMs are not complete, or circumstances where only a few QNMs 
are known and others are not. (This is in practice often the case, since QNMs with large 
— Im uj are difficult to obtain, even numerically. The case of gravitational waves propagating 
away from a black hole is precisely of this type.) Thirdly, QNMs, being complex, are not in 
general plagued by nodes for real x, where f(x) would be singular. In the case of NMs, the 
zeros of excited state wavefunctions require special treatment fl21~| , |23ll ; fortunately for QNMs, 
this is not a problem. 



In terms of f(x), ( |1.3|) becomes the Riccati equation 
— rice 



f'(x) + f{x) - [V (x) + fiV^x)] + uj 2 = . (2.1) 

The function f(x) satisfies both the left boundary condition (f(x) — > —iuj as x — > —00) 
and the right boundary condition (f(x) — > +iu as x — > +00), with uj in ( |2.1| ) being the 
eigenvalue uj. However, for any uj (whether or not an eigenvalue), we can define two solutions 
f±(uj,x) by the boundary conditions f±(u,x) — > ±iu as x — ► ±00. At an eigenvalue uj, 
f+(uj, x) = f_(u,x) = f(x). We expand the eigenvalue uj and the eigenfunction / in powers 
of \i: 

— omega 
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UJ = LU + fJXVl + /i 2 Co> 2 + 



(2.2) 



/ = fo + 9 = fo + H9\ + ^ 2 92 H , (2.3) 

where / , assumed known, satisfies the Riccati equation fl2.1|) with the potential V and 
frequency loq. 

Putting (|2.2|) and (|2.3|) into the Riccati equation (|2.1|) and upon comparing powers of /i, 
one finds, in a straightforward manner, that 

— gn 

g' n + 2f g n + 2uj uj n = V n (2.4) 

for n = 1,2, in which V\ is the perturbing potential in ( |1.3| ), and V n , n > 1, is the 
effective nth order potential, defined in terms of a combination of lower-order quantities: 

— vn 

n-l 

V n {x) = [9i{x)g n -i{x) + COiUJn-i] ■ (2.5) 

1=1 

Using the integrating factor exp[2 / dxfo(x)}, ( |2.4D can be solved; but since the resultant g n , 
related to an eigenfunction, must satisfy two boundary conditions, this imposes a condition 
on the remaining free parameter uj n . The details of the derivation will be given elsewhere 



24 1 . Ignoring convergence problems for the moment, and noting that the logarithmic deriva- 



tives at spatial infinity are unaffected by the perturbation, by definition of the boundary 
conditions, one obtains the following formal expressions, which constitute the core of LPT 
2|: 
-result 

2u; o (0o|0o) 

Here we have introduced the suggestive notation 
— matelm 

/•CO 

(0o|K|0o) = / V n (x)$(x)dx (2.7) 



-norm 



oo 

2 



(0 o |0o) = / P {x)dx . (2.8) 

These results express the nth-order correction to the eigenvalue in quadrature in terms 
of lower-order quantities. The nth-order correction to the logarithmic derivative is given by 
— function 

g n {x) = ( f dy [V n (y) - 2u u n ) <t> 2 {y)) <j>f{x) . (2.9) 
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When g n is substituted back in (|2.5|) to obtain V n+ i, the cycle of iteration is complete, 
and one can in principle obtain the corrections to any order. The strategy outlined here is 
basically the same as in the conventional LPT for hermitian systems. 

The above expressions are formal and suggestive, but hide an essential problem: because 
the unperturbed wavefunctions go as e ° x as x — > ±00, the integrands contain factors 
(j)o(x) 2 oc e 270 '^ as Id 



00 (70 = — Im ujq), rendering the integrals ( p.7|) and (|2.8|) divergent. 
It is therefore critical to handle the asymptotic regions carefully, which, as we shall see, 
may be regarded as a way of giving meaning to or regularizing these formally divergent 
expressions. 

One can deal with the asymptotic regions in several slightly different ways ||24|| , and we 
here give a simple version that relates most directly to the actual numerical implementation 
of the program. To do so, we return to ( |2.4| ), insert the integrating factor exp[2 J dxfo(x)] = 
<f>o(x), and integrate from L_ to L + : 

— intl 



4>l(x)g n (x] 



l + 2u uj n 



dx 4>l(x) 



dx 4>l(x) V(x) 



(2.10) 



(Had we taken L± = ±00 instead of finite values, we would return to the formal result stated 
above.) 

Now it has to be noticed that g n contains implicitly the unknown nth order frequency 
shift uj n . To see this, we note that g(L±) is to be matched to f±(w, L±)—f (L±); consequently 
its nth-order part is 

— gnmatch 



V n 9n = [/±(w,L±) - fo(L±)} n = [f±(u,L±) - f±{uJ ,L ± )} 



f±(uQ + fJWx-\ hfi n u n -i,L±) - f±(u ,L±) 



d 

+ fi n uj n — f±{uJo, L±) 

n OUJ 



d 

A± n + u n —f±(u , L±) 



(2.11) 



where the subscript on the square brackets indicates taking the nth-order part, and we have 
separated out the term that depends on oj n . The quantities A± n express the difference 
between two logarithmic derivatives at the asymptotic points. Thus, collecting the terms 



that do and do not contain u> n , we obtain from ( |2.10| ) again an expression like (|2.6|), except 
that the matrix element and norm are now given by 
— matelm2 



(0O|K|0 O ) 



V n {x)(t)l{x)dx + A_ n 2 (L_) - A +n (f) 2 (L + ) 



(2.12) 



— norm2 
(4>o\4>o) = 



L+ I 

4>l(x)dx + — 



d 

%(L+)-^f+(uj ,L A 



(2.13) 



Thus, provided we have a way of calculating the values of f±(uj,L±) at the asymptotic 
positions L±, these alternate expressions ( P-12j ) and (p.13 ) then provide practical ways of 



7 



evaluating the matrix element and the norm, without infinite integrals or divergent expres- 
sions. In other words, the difficulties alluded to earlier have been eliminated. The choices of 
L± are arbitrary, and the shifts uo n must be independent of these choices; in fact, it is easy 
to see that ( ]2.12|) and (|2.13|) must be separately independent of L±, simply because (|2.13[) 
relates only to the unperturbed state, whereas ( |2.12| ) depends on the arbitrary perturbation. 
This independence can be verified both analytically by differentiating the expressions with 
respect to L± and using the Riccati equation, or numerically by evaluating the expressions 
for different values of L±. 



III. GENERAL PROPERTIES 

The properties of the perturbed QNM's can be discussed at three levels of generality: 
(a) properties about open systems in general (in contrast to conservative systems); (b) 
properties for any perturbation of a Schwarzschild black hole; and (c) results for a specific 
perturbation. This Section deals with the first two, and a specific model of perturbation is 



presented in Sect. IV. Although the LPT is valid to all orders, for the present application 



we focus on the first-order shift. 



A. Open systems in general 

The result in ( |2.6|) has been written in a way formally similar to the hermitian case. 
The factor 2u Q occurs because the eigenvalue is u 2 rather than u. Since the numerator 
and the denominator in ( |2.6|) are separately independent of L±, they can be given physical 
interpretations as a generalized matrix element and a generalized norm respectively. 

The generalized norm has been introduced previously in a restricted form |26[ applicable 



only to cases where the potential has "no tail"; in that case A± can be obtained readily, 
since f±(u>,L±) = In that restricted form, the perturbation theory for the QNMs of 

the wave equation has been developed and applied [^7] . For systems such as black holes, the 



"no tail" condition must be removed, as sketched in Sect. II. In either case, the generalized 
norm has some unusual properties; (a) It involves 0q rather than |0o| 2 , and is in general 
complex, (b) It involves surface terms at x = L±, though the value of the entire expression 
is independent of the choice of L±. Thus, it is not a norm in the strict sense, but rather a 
useful bilinear map. In cases where the system parameters can be tuned so that the leakage 
of the wavefunction approaches zero (e.g., Vq(x) contains a tall barrier on both sides), then 
the generalized norm reduces to the usual (real and positive-definite) norm for a NM. 

It is useful to define a function H(x) for each QNM which depends only on the original 
unperturbed system: 

-hdef 

5U ~ H{x) = M *l : . (3.1) 



Both the magnitude and the phase of H(x) are well defined and physically significant. The 
magnitude implies that we can now give a precise meaning to the normalization of a QNM, 
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even though the wavefunction diverges at infinity. The phase of H(x) determines the phase 
of the first-order shift u)\ for a real and positive localized perturbation V\{x). The phase is 
intriguing because it has no counterpart for a hermitian system — in that case, H(x) must 
be real and non-negative. For an open system, H(x) can be negative or indeed have any 
phase, so that a positive perturbation could lead to a decrease of the frequency. 

The functions H(x) are then convenient objects for discussing the effect of any pertur- 
bation on the QNMs of a given system. We next present some properties of H(x) for the 
Schwarzschild black hole. 



B. Schwarzschild black holes 

First consider the normalization of these QNMs, made possible by the introduction of 
the generalized norm. Since ( |2.8j ) is independent of L±, it is particularly convenient if we 
take L_ = L + = 0: 

-defA 

(0ol0o) = 2^; 0o(o)2 {L f+{uJOi0) ~ L f - {uJOi0) } = ° (o)ir • (3,2) 

The logarithmic derivatives at x = can be obtained by integrating the differential equations 
from x = ±00, and are readily calculated making use of Leaver's solution |2£| (see Appendix 
|B|). The parameter K is a convenient way of expressing the normalization, and Table | lists 
the values of Kj (2M) for the lowest few QNMs (labeled by j) of the Regge- Wheeler potential 
for each angular momentum I. They increase as j increases and show clear patterns, e.g., for 
I = 0, s = 0; K is large for odd j, while for I = 2, s = 2; K is large for even j. It is interesting 
to ask in what way the values of K can characterize the Regge- Wheeler potentials, and to 
investigate their meaning beyond the perturbation calculation (noting that K has nothing 
to do with the perturbation). 

Next consider the functions H(x). Fig. ^| shows Re H(x) and Im H(x) for I = 1 scalar 
waves, versus the coordinate x/(2M). The diagrams refer to the lowest QNMs (labeled as 
j = 0, 1, ■ • ■ , 5) of that angular momentum. Since the wavefunctions grow exponentially as 
\x\ — > 00, we have multiplied H(x) by exp[— 2T^Jl + (x/2M) 2 ] (r = — Im 2Mujq) for ease of 
plotting. These figures show an essentially plane wave behavior for large \x\, but the phases 
are nontrivial. We note that for a localized perturbation /jlVi(x) = fi5(x — Xx), the first- 
order shift in frequency u>i is given by H(xi), and therefore can be read out directly from the 
figures. For example, Table |T| lists the magnitude and phases of H(x) for I — 1, s — and 
the lowest few j's at a fixed value x = x\ (r = ri), for x%/2M = 0, 5, 10 (r\/2M = 1.28, 
3.93, 8.05); r x /2M = 1.02 and 1.08 {x x /2M = -2.89, -1.45). The patterns are different for 
these values of x% and not simple, demonstrating that a localized perturbation will push the 
QNMs along different directions in the complex frequency plane. 

There is a more complicated structure in H(x) for smaller x, best exhibited if the same 
data are plotted versus r/2M, as in Fig. [3[ (The region r/2M < 2.0 corresponds to x/2M < 
2.0.) The many oscillations for x < are compressed into a small region near r/2M = 1 
and are not visible in this plot. One interesting feature is that both Re H(x) and Im H(x) 
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alternate in sign as (— 1) J near the event horizon. The magnitude and phase of H(x) at 
the sample positions r = n where rj2M = 1.02 (x 1 /2M pa -2.89) and n/2M = 1.08 
(x\/2M pa —1.45) are also shown in Table 

The results here described depend only on H(x), i.e., on the properties of the unperturbed 
system. These examples indicate that localized perturbations can generate a rich pattern 
of frequency shifts (in contrast to shifts all of the same phase in the case of the NMs of a 
conservative system). In turn, this means that there are much better prospects of learning 
something about the perturbing potential from the observed shifts. 

Of course, the richness of the pattern could be diluted if the perturbation is not localized 
in x, but has a spatial extent Ax large compared to the typical wavelength A of oscillation 
of H(x), A pa 27r/Re ujq ~ a few M (see Fig. ||]for example). In the next Section we evaluate 
and discuss the perturbation that arises from the presence of a dust shell, which leads to a 
perturbation potential extending over a range of x. 



We consider a Schwarzschild black hole surrounded by a static spherical shell of matter 
at a certain radius outside the hole. The equation of state of the matter making up the shell 
is chosen so that the shell remains static at the given radius. This example demonstrates 
that even such a simple model can lead to intriguing features in the QNM spectrum of the 
perturbed system. 

There are two mass parameters in the black hole plus shell system: the total ADM mass 
of the system measured at infinity M Q , and the mass of the black hole as measured by its 
horizon surface area M a (i.e., the surface area of the event horizon is 16tcM 2 ). In the limit 
of M a — > M a , we return to the unperturbed case of a bare black hole. Thus the parameter 
\i = (M — M a )/M a is a measure of the perturbation. The perturbation further depends on 
the circumferential radius r = r s where the shell is placed. We shall study the shifts as a 
function of \i and r s . 

A static and spherically symmetric metric can be written as — metric 



IV. A MODEL PROBLEM 



A. Description of the model 



ds 2 = -A{r)dt 2 + B(r)dr 2 + r 2 {d6 2 + sin 2 6 2 dip 2 ) . 



(4.1) 



As usual, we define the mass function m(r) such that 



B(r) = 1 



2m(r) 



(4.2) 



Qrr 



r 



The functions m(r) and A(r) satisfy the equations 



dm 
dr 







r ^ r s 



(4.3) 
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2A dr r(r — 2m 
Hence we have 



1 dA m 

(4.4) 



f \ J M ° r > T » (A K\ 

mir) = < (4.5) 
w I M a r < r s y ' 



1 - 



2M Q 



r > 7% 



A(r) 



a 1 



2M n 



(4.6) 



r < r. 



where 



a 



1 - 2M Q /r s 
1 - 2M a /r s 



(4.7) 



The metric ( |4.1|) becomes — metric2 



ds z 



2M n 



-a 1 



2M a 



rft 2 + 1 - 



2M a 



rfr 2 + r 2 (d6 2 + sin 2 



<ir 2 + r 2 (d6 2 + sin 2 



r > 7% 



r < 7% 



(4.8) 



The physical meanings of the parameters r, M D and M a are now clear from ( |4.8| ). We see 
explicitly that in the limit M Q — >■ M a , the value of a reduces to unity, and Q4.8Q reduces to 
the familiar Schwarzschild metric. 

Next we consider waves propagating on this perturbed Schwarzschild background. In 
this paper for simplicity we focus on the case of scalar wave (s = in Q1.2Q ), described by 

d,d^ = [y/^ggT^) = • (4.9) 

To obtain the Klein-Gordon equation suitable for the perturbation analysis, we separate ip 
into the radial and angular parts: 



oo / 

V>M) = E E Ri{r,t)Y lm {e,cp) . 

1=0 m=-l 



(4.10) 



The radial function satisfies 



2A A' AB'\ n/ Al(l + 1 



B \rB 2B 2B 2 



R'i- 



-Ri = d 2 Ri 



(4.11) 



where ' = d/dr. Next introduce the tortoise coordinate 

^B 



x 



dr 



(4.12) 
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to push the event horizon (r = 2M a ) to — oo. We leave the zero point of x to be specified 
later. The equation for Ri becomes 



d^Ri H — \l — d x Ri — — -Ri = d^Ri 



r V B 



(4.13) 



To cast the equation into the standard Klein-Gordon form, we introduce a function 
(hereafter the label I will be suppressed) such that 



R(x,t) = h(x)$(x,t) 



(4.14) 



with 



d x h + -J^h = . 
r V B 



(4.15) 



This function $(x,t) then satisfies the Klein-Gordon equation ( |1 . 1|) with the effective po- 
tential 



T _, , Al(l + 1) A' AB' 
V(r) = — H 

Substituting in the expressions for A and B, we get 



(4.16) 



V(r) 



1 



2M n 



/t r (5(r -r s ) 



1(1 + 1) 2M 



+ 



= V sc (M ,r) 



a 1 



2M n 



Z(Z + 1) 2M a 

9 T Q 



al/ sc (M a ,r) 



r > r s 
r = r s 
r < 7% 



(4.17) 



where Vsc(A^, r) is the standard unperturbed potential for a Schwarzschild black hole with 
mass M (i.e., equation (|1.2p with s = 0). The strength of the 5-function at the position of 
the shell r s is given by 



lim 



V(r) dr 



1 



2M \ [iM a 



(4.18) 



The remaining degree of freedom in setting up the perturbed problem is the zero-point 
of the tortoise coordinate x. We choose it so that the effective potential is the standard 
unperturbed form outside the shell, in terms of the ADM mass. The zero point is then 
determined by the requirement of continuity of x. 



x 



r + 2M ln(r/2M - 1) 

1 [r + 2M a ln{r/2M a -l)]+x c 
'a 



r > r s 
r < r. 



(4.19) 



where x c is the constant which makes x continuous. 

Thus we arrive at the perturbation problem ( |1.3j ) as discussed in the sections above with 
x being the tortoise coordinate, 
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V (x) = V sc (M o ,x) (4.20) 



, \ J for a; > x s , . 

= I ^ _ x j + a ^ 5C (M a , r) - V^c(M , r) otherwise ' ^ ' 



where 



and 



■>■, =U 2.U ln(^--l) (1.22) 



« = Hm / cfe = — / 4 _23) 



This perturbation consists of two parts: a 5-function at the position of the shell, plus a 
contribution inside the shell extending all the way to the horizon (x — > — oo or r — > 2M a ). 
There is no perturbation outside the shell. 

We should make one further comment concerning the setup of the problem. As the 
equation of state of the matter shell is chosen to make the shell static, one need to ensure 
that the equation of state satisfies the energy conditions. It is straightforward to show that 
the dominant energy condition is the first to be violated if the shell is placed too close to the 
horizon, as shown in Appendix [A|. Thus, the discussion will be limited to the perturbations 
that do not violate the dominant energy condition. 



B. Numerical solution 



The full potential V can be cast into a standard Regge- Wheeler potential for both x < x s 
and x > x s , and hence exact QNMs can be obtained by numerically evaluating the Leaver's 
solutions of the Regge- Wheeler equation (see and Appendix |TJ). The numerical solution 
will be used to examine the accuracy and validity of the perturbative result. Details of the 
numerical scheme are sketched as follows. 

Firstly, as there is a 5-function in the effective potential V(x), the logarithmic derivatives 
at x = x s are related by 



lim 



<)'(u,x 8 + e) 4>'(u,x s -e) 



(u>,x s + e 



(uj,x s 



(4.24) 



Secondly, outside the shell, the effective potential is the Regge- Wheeler potential with 
mass parameter M and the QNM wave function satisfies the outgoing wave boundary 
condition at the spatial infinity x — * oo. Hence, 



where 



<f)(uj, x) = 0+(2M o w, x/2M ) , 
(u, x) is the outgoing wave solution of the scaled Klein-Gordon equation 



4 + ^ - V(r) 
ax z 







(4.25) 



(4.26) 
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with 



and 



V(f) 



1(1 + 1) 1 



(4.27) 



x = f + ln(f - 1) . (4.28) 

Note that the quantities with tilde are dimensionless. 

Thirdly, inside the shell, it is also easy to show that the wave function satisfies 

2 



+ 



OJ 

fa 



-V sc (M aj r) 



x < x* 



and 



Hence 



x = r + 2M a In 



2M n 



\fa(x — x c 



{uj,x) = </>.(2M a u/y/a,x/2M a ) . 



(4.29) 



(4.30) 



(4.31) 



Here <p-(u}, x) is the ingoing wave solution of ( f4.26|) at x, — > —oo. 

Therefore, by connecting the two logarithmic derivatives at x = x s with (|4.24j ), the QNM 
condition can be written as 



(f)' + (2M u,x) 



(f) + (2M uj,x) 
where p = -faM jM a , 



p<p'_{2M uj/p,x) 
i= Is / 2 M„ <p-{2M u/piX) 



2M k 



x=x 3 /2M a 



x s = r s + 2M a In 



\2M a 



i^j = y/a(x s - x c ) 



(4.32) 



(4.33) 



and the prime represents the derivative with respect to x. The functions <ft + and 0_ can 
be computed by the Leaver's solutions of the Regge- Wheeler equation (see [28| and Ap- 
pendix 0). Therefore, exact QNMs can be obtained by solving the nonlinear equation ( [4.32| ) 
by standard root-searching methods. 



C. Dependence on fi and convergence 



Given the perturbation function fiV\ in (|4.21| ), the shifts u\ and u 2 can be evaluated by 



( p.6| ). The detailed treatments will be given in Appendix |C|. 

Fig. ^ shows the 0th, 1st and 2nd order perturbation results for / = 1 scalar waves 
together with the exact numerical results for a black hole plus shell system. The parameters 
are [i = 0.02, r s = 2.52M a (x s ~ 0). This figure shows that the perturbation formalism does 
give the correct shifts for small fi. 

To see the convergence more clearly, we plot in Fig. |5| the magnitude of the error in the 
frequencies in the 0th, 1st and 2nd order results versus [i. The plot shows the case of / = 1, 
s = 0, j — 1 (first excited state), with the shell located at r s = 2.52M a . The error of the 
nth-order result goes as u n+1 , as it should. 
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D. Dependence on shell position 



We next study the dependence on the parameters of the shell. Fig. |^ shows the trajec- 
tories of the lowest damping QNMs (j = 0, 1, ■ • • 6) as the position of the shell moves away 
from the event horizon. The plot shows the case of I = 1 scalar waves with fi = 0.01 based 
on numerical calculation. The regions near the j = and j = 1 modes are shown in greater 
detail in the insets. The calculation is terminated when the QNM is near the imaginary uo 
axis, since it is difficult to perform the numerical calculation with sufficient accuracy in that 
region. It is seen that there are interesting changes in signs with increasing j. Some QNMs 
move toward the imaginary u axis, and the higher damping QNMs either move toward the 
imaginary uj axis or move upward toward the origin. It is also seen that when the shell 
is placed far away from the event horizon, the QNMs move away from their unperturbed 
positions, with the higher damping modes moving with higher speed. This behavior can be 
understood from the perturbation formula (|2.6|). When x s /2M a 3> 1, we have, according to 



(fO|) and (TO , — asyom 

u x ~ e 2luJ ° Xs /x 2 for x s /2M a > 1 . (4.34) 

Hence QNMs move more rapidly away from the unperturbed positions when the shell 
is far away from the event horizon (x s /2M a ^> 1) and the higher damping modes 
(— Im 2M a ujQ ^> 1) move with higher speed. 

Fig. [7] shows the first-order perturbation result for the j = mode for different values 
of x s . At large x s , the trajectory shows a spiral structure, which can be explained from the 
first-order perturbation formula. According to ( |4.34j ), 

du dui e 2iuJ ° x ° . 1 , ao . 

— « -j^ — for x s > 1 . (4.35) 

The exponential factor e 2luJ ° Xs leads to the spiral structure. For larger j modes, there is no 
such spiral structure (Fig. |6]) because the higher-order corrections become large before the 
spiral appears, as shown in Fig. |j. 



V. CONCLUSION 

We have applied the LPT for the QNMs of open systems to the study of gravitational 
waves propagating away from black holes. QNM gravitational wave signals from black 
holes will be detected soon, and many black holes are expected to be perturbed by their 
astrophysical environment, e.g., by an accretion disk. The study in this paper on scalar 
waves represents the first step in the study of the QNM spectra of waves propagating in 
a dirty black hole background, perturbed by the astrophysical environment. Its further 
development can be of interest to gravitational wave astronomy, among other applications. 

Although the QNMs of any system can in principle be obtained through brute force nu- 
merical integration, it is nevertheless highly desirable to have a convenient set of perturbation 
formulae that one can understand the system with. We note that the Rayleigh-Schrodinger 
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perturbation theory for NMs is tremendously valuable, even though the NM systems are 
more easily handled numerically than QNM systems. 

We have shown in a simple example that a perturbed black hole spectrum can have 
interesting features, and demonstrated how these features can be understood with the per- 
turbation formula. In summary, we raise the importance of studying QNMs of dirty black 
holes, and show how it can be done in a perturbative formulation. 
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APPENDIX A: DOMINANT ENERGY CONDITION 

The shell cannot be placed too near to the black hole; otherwise, the dominant energy 
condition 

1 5/| > \Sg 4 1 (no sum on q) (Al) 
will be violated, where q denotes any unit specelike vector, and 

S/= hmj^ T^df . (A2) 

It is straightforward to compute all the non-trivial components of the Einstein tensor for 



the metric (fi~T 



A' 1 1 

g/ = 7ab + ^b~^ (A3) 

c e = cv = J± A' 2 B' A'B' A" 

9 v " 2rAB 4A 2 B 2rB 2 AAB 2 2AB 1 ' 

c t = 1 1 B ' ( A5 ) 

,y2 ty 2 ty 2 

where ' = d/dr. The component is non-zero only if it contains terms involving the 
second derivative of A or the first derivative of B. Hence 

S { * = = — lim / — — dr (A6) 

* 87T e^0+ Jr s -e tB^ 2 V ' 



2nr 2 yv , s v 
using the Einstein equations = 87tT^. Similarly, we have 
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1 - M Q /r s 1 - M a /r s 



57rr. 



1 - 2M /r s Jl - 2M a /r 



(A8) 



[I]) together with (|A7|) and (|A8|) imply a lower bound on r s for a given p. For /i = 0.01, 
it is found that the minimum r s is 2.26 M a . 



APPENDIX B: LEAVER'S SOLUTIONS 



In the analytic study of the Regge- Wheeler equation, it is convenient to rewrite the 
Regge- Wheeler equation in terms of the coordinate r: 



r(r — 11 



+ </>.r - 



£i + j( J+ i)+£ 

r — 1 r 







^Bl) 



where p = —iu, £ = 1 — s 2 . For simplicity, we have set 2M = 1. The equation flB~T| ) can be 
transformed to the generalized spheroidal wave equation [p8fl . For an arbitrary frequency tu, 
we define two solutions (f>±(uj,r) which satisfy the boundary conditions 



<j>± oc e ±iuJX for x -> ±oo 
The analytic solutions of these two functions are given by p8 



n=0 

oo 



(u;, r) = r 1+s (r 



- iy e - pr J2 a «(2p + l)n + 1 + 2p + n, 2s + 1, 2pr) 



n=0 



where 



(2p + 1)„ = 



T(2p + l + n) 

r(2 P + i) 



(B2) 

(B3) 
(B4) 

(B5) 



is the Pochhammer's symbol, U is the irregular confluent hypergeometric function [p0|] , and 
a n is determined by the following three-term recursion relation: 



with 



a n a n+1 + (3 n a n + 7 n a n _i = n = 1, 2, 



a n = for n < 



(B6) 



(B7) 



The value of ao is arbitrary and is related to the normalization of <p±. The quantities a n , 
j3 n and 7 n are defined to be 



a n = (n + l)(n + 2p+l) 

/3 n = - [2n 2 + (8p + 2)n + 8p 2 + 4p + 1(1 + 1) + ^ 
7 n = n 2 + 4pn + 4p 2 + ^ - 1 . 



(B8a) 
(B8b) 
(B8c) 
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In most cases, we only need the logarithmic derivatives 

1 d(j>± 



f±(v,x) 



dx 



(B9) 



While f -(id, x) can be easily obtained from (p3|), the calculation of f+(uj, x) is less straight- 
forward as it involves a sum of irregular confluent hypergeometric functions, which are 



notoriously different to evaluate |28| . In the following, we develop a numerical algorithm to 
calculate f + (ui,x). 

We define two sequences 



= (2p + l) n U(s + 1 + 2p + n, 2s + 1, 2pr) 

T« = -^±l±g- M( 8 + 1 + 2p + n, 2s + 1, 2pr) 
r(2p + n — s + 1) 

where M is the regular confluent hypergeometric function. Hence we have 

<f>+(w,r) = r 1+s (r - Ife-^h^u, r) , 

where 

oo 

r ) = H a nTt ] ■ 
n=0 



(BIO) 
(Bll) 

(B12) 

(B13) 



f+(u,x) can then be expressed in terms of h as follows: 

"1 + s P 



1 



r J <p + dr 



i-I 

r 



+ 1 d/i 

r — 1 ^ ft, c?r 



(B14) 



It can be shown, from the recursion formulae of the confluent hypergeometric functions 
330, that 



(2p + n)(2p + n + l)r£l - (2p + n + 1)(1 + 4p + 2n + 2pr)T^ 
+ (s + 1 + 2p + n)(2p + n - s + 1)T^\ = 



and 



dJf^ _ s + l + 2p + n 
dr r 



2p + n + l-s ( _) _ ( _ } 



2p + ?2 + 1 



(B15) 



(B16) 



For large n, T^> approaches the asymptotic expressions [ pi 
v^Fr(2p + l + n) 



- — 

T(2p + 1) T(2p + 1 + n + s) \2pr) 



(2npr)- 1/4 e pr exp 



-2(2npr) 1/2 



2^7rr(2p+ 1 + s + n) \2pr 



(2npr)~ 1/i e pr exp +2(2npr) 



,1/2 



(B17) 



(B18) 



The general solution of the difference equation ( |B15|) is a linear combination of It 
follows from ( |B17| ) and ( [B18[ ) that T^ + ) increases with n, while decreases with n. 

Hence any linear combination of will eventually be dominated by the term containing 



T^ + ) as n — > oo. The computation of T„ +1 from ) and T„_{ is therefore unstable. The 
recursion relation ( |B15| ) must be used in the reverse direction. In practice, we choose an 



arbitrary value of Tjy ^ (N ^> 1) and set 



T 



(-) 



T 



(-) 



N-l 



N-l 



(-) 



xT 



(-) 



N 



■ N 



^ } 2p + N + s fN-l\ s - 1/4 



N 



2p + N 



N 



exp 



2(2 pr) 



1/2 



N + y/N-l_ 



(B19) 



The values of for n < N — 1 are then calculated by the recursion formula ( |B15p . 
determined in this way will differ from its original definition in ( |B10| ) by a multiplicative 
constant. This does not concern us since we are only interested in the logarithmic derivative 
of + , which is independent of the constant factor and is readily evaluated by ( |B 14|) . 

Leaver's solutions can be applied to compute the logarithmic derivatives fo(L±) in (|2.13|) 
and to evaluate the QNMs of the shell model of dirty black hole described in Sect. [TV]. 



APPENDIX C: EVALUATION OF THE PERTURBATION FORMULAE 



For simplicity we focus on the first order perturbation calculation, which is given by 

(0 O |Wo) 



2u> o (0o|0o) 



(CI) 



We have discussed the evaluation of the generalized norm (0o|0o) i n Sect. [IIIB| . It suffices to 
evaluate the generalized matrix element (0o|14|0o)- F° r Vi given by Q4.21|) , the generalized 
matrix element is 



<0o|Wo> 



Vi(x)4>l(x) dx + K(f) (x s ) 2 + 0q(L_)— f-(u , L. 



(C2) 



Ai=0 



where L_ is any real number smaller than x s . It has been argued earlier that this expression 
is independent of the choice of L_, and for any finite L_ the expression is finite. Thus the 
entire perturbation scheme requires no regularization. 

However, the last term in (p2|) involves the change of /_ due to the presence of V\. 
This calculation can be bypassed by the following trick. First, since the entire expression is 
independent of L_, we push L_ to — oo and define a function 



Z(u) 



where 



lim 

L— — oo 



d 

(w,r*)yi(r»)dr, + 4> 2 _(u,L_) — f_(u,L_ 



ti=0 



+ K(j)l_(uJ,X s 



(C3) 



is a function satisfying the unperturbed Regge- Wheeler equation with the 
ingoing wave boundary condition on the event horizon. This expression as a whole is finite 
for every uj, but for lm.(u) > 0, each term is separately finite; in fact the surface term 
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vanishes as L_ — > — oo because the factor 0q_(u;,L_) oc e 2lujL - now decays exponentially. 
Evaluating in this domain, we have 



' 2 



Z{lo) = I 00-^) r *)^i( r *) dr* + K0 o „(tu,x s ) for Im(u;) >0 . (C4) 

J — oo 

The prescription is therefore to evaluate in this domain and analytically continue to ujq: 

(0o|Vi|0o) = Z(u ) = / (^{r^V^r^dr* + K(j>l(x s ) , (C5) 



oo 



We note that this merely provides an alternate evaluation of an expression that was mani- 
festly finite to start with. 

The integral can be easily carried out analytically by noticing that the potential 
Vsc(M, x) in ( [4.1 7| ) can be expressed as a sum of exponentials at negative x; the calcu- 
lation is sketched below and the details of the straightforward arithmetic can be found in 
Ref. Q: 

^(*.*> = dp v "G-2Sr) (C6) 

OO 

V sc (l/2, x) = J2 Cke ki for ~x < , (C7) 



k=l 



where 

c k = 1(1 + l) 73 , fe + (1 - s 2 ) litk (C8) 

m 

7n,m = e~ m p(p, n, m) . (C9) 

p=0 

Here /i(0, n, 1) = 1 — n, n, 1) = n, [i(p, n, m) = for p < and p > m, and 

in + p + m — 2) u(p — 1, n, m — 1) + (m — 1) uOp, n, m — 1) 

u(p, n, m) = for m > 2 . 

m 

(CIO) 



It can be shown that the unperturbed wave function O can be expressed as 29 



0o(x) = £4Me fe * /{2Mo) , (Cll) 

k=0 

where 

i fc-i 

d fe( w o) = 771 7T- — V dk(wa)Ck-m , (C12) 

k(k - 2iuj ) ^ 

and do is a constant related to the normalization of the wavefunction. 

It is obvious that V\ given by ( 4.21|) can also be expressed as a sum of exponentials 
plus a 5-function. Hence the integration in (|U2[) can be performed analytically, giving the 
generalized matrix element and the first-order shift uj\. 
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FIG. 1. The distribution of QNMs of a Schwarzschild black hole for / = s = 2. 
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FIG. 3. Graph of Re H{x)e' 2T ^ l+ ^ x / 2M ^ (solid line) and Im H(x)e~ 2r V 1 +( x /2M) 2 (d as h line) 
vs r/2M. 
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FIG. 4. The perturbed QNM spectrum of a black hole plus shell system for I = 1 scalar 
wave with fi = 0.02 and r s = 2.52M a . The horizontal axis is Re 2M q uj and the vertical axis is 
Im 2M„w. The 0th, 1st and 2nd order perturbation result is indicated by circles, squares and 
triangles respectively. The exact numerical result is represented by stars. The diagram on the 
upper left shows the distribution of QNMs of the lowest damping modes (from j = to j = 6). 
The other diagrams are the magnification of the region around each mode. 
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FIG. 6. The trajectory of the lowest QNMs of I = 1 scalar waves for /i = 0.01 and r s /M a 
varies from 2.26 to 60 based on exact numerical calculation. The circles are QNMs of a bare 
Schwarzschild black hole with mass M Q ; the squares are the QNMs for r s = 2.2QM a (The dominant 
energy condition is violated when r s < 2.26M a ); the triangles show the positions of QNMs at 
r s /M a from 6 to 60 in intervals of 6. The upper and lower inset show the regions near j = and 
j = 1 in details respectively. We stopped the calculation when the QNMs approach the imaginary 
lo axis because it is difficult to compute it accurately in that region. Thus, some triangles appear 
"missing" for some higher damping QNMs because they have already moved to the imaginary lo 
axis. 
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FIG. 7. Trajectory of j = mode based on exact (solid line) and first order (dash line) cal- 
culation. The zeroth order result is indicated by a circle and the square represents the mode for 
r s = 2.26M a . The triangles show the mode for the positions of the shell at r s /M a from 6 to 36 in 
intervals of 6 for exact result, and the plus symbols show the first order result for the shell at the 
corresponding positions. 
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FIG. 8. Same as Fig. [7| but for j = 1, 2, • • • 6. For these j, the perturbation breaks down before 
a spirl forms. 
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TABLES 



TABLE I. The normalizing factor Kj (2M) (expressed as magnitude and phase) for the lowest 
QNMs of the Schwarzschild black hole with mass M. 





I = 0,s ~- 


= 


I = l,s- 


= 


l = 2,s 


= 2 


J 


magnitude 


phase 


magnitude 


phase 


magnitude 


phase 





2.4461 


36.5° 


4.5155 


36.7° 


5.5397 


18.7° 


1 


193.24 


18.4° 


20.544 


68.1° 


6.9371 


82.7° 


2 


6.9935 


12.9° 


12.319 


22.7° 


22.770 


-66.9° 


3 


894.17 


-10.2° 


33.553 


10.1° 


9.4942 


39.4° 


4 


12.059 


8.4° 


18.945 


20.5° 


63.958 


-89.7° 


5 


2279.0 


-10.4° 


60.316 


-14.2° 


14.358 


25.8° 


6 


17.087 


6.7° 


25.087 


26.6° 


66.779 


-100.9° 
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TABLE II. The value of H(x) (expressed as magnitude and phase) for several sample positions 
at x = xi (r = ri), for I = 1 scalar waves propagating on a Schwarzschild black hole. 





x x /2M = 


= 


x ± /2M - 


= 5 


Xl /2M -- 


- 10 


n/2M = 


1.02 


n/2M = 


1.08 


3 


magnitude 


phase 


magnitude 


phase 


magnitude 


phase 


magnitude 


phase 


magnitude 


phase 





2.21 x TO" 1 


-36.7° 


9.18 x 10- 1 


-122.2° 


6.44 x 10° 


-147.1° 


4.16 x 10- 1 


4.8° 


2.66 x 10- 1 


-38.7° 


1 


4.87 x 1(T 2 


-68.1° 


5.23 x 10 1 


109.8° 


2.40 x 10 4 


52.1° 


2.69 x 10 1 


-141.1° 


2.14 x lO" 1 


150.9° 


2 


8.12 x 1(T 2 


-22.7° 


2.88 x 10 3 


-68.0° 


1.28 x 10 s 


-165.3° 


3.14 x 10 1 


39.1° 


1.03 x 10° 


-28.5° 


3 


2.98 x 1(T 2 


-10.1° 


2.27 x 10 5 


96.8° 


1.60 x 10 12 


-30.4° 


4.20 x 10 2 


-139.9° 


3.32 x 10° 


156.4° 


4 


5.28 x 1(T 2 


-20.5° 


2.28 x 10 7 


-99.7° 


2.51 x 10 16 


112.4° 


6.00 x 10 3 


40.2° 


1.19 x 10 1 


-19.2° 


5 


1.66 x 1(T 2 


14.2° 


2.61 x 10 9 


65.9° 


4.56 x 10 20 


-97.0° 


9.01 x 10 4 


-140.2° 


4.27 x 10 1 


164.0° 


6 


3.99 x 10" 2 


-18.0° 


3.22 x 10 11 


-126.3° 


8.91 x 10 24 


59.6° 


1.40 x 10 6 


39.4° 


1.57 x 10 2 


13.6° 
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